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C^l ■ Abstract. The Hamiltonian for dynamic geometry generates 

the evolution of a spatial region along a vector field. It 
includes a boundary term which determines both the value 
of the Hamiltonian and the boundary conditions. The value 
1 gives the quasi-local quantities: energy-momentum, angular- 

. momentum/center-of-mass. The boundary term depends not 

only on the dynamical variables but also on their reference 
values; the latter determine the ground state (having vanishing 
quasi-local quantities). For our preferred boundary term 
for Einstein's GR we propose 4D isometric matching and 
extremizing the energy to determine the reference metric and 
connection values. 



1. Introduction 

Energy-momentum is the source of gravity. Gravitating bodies can 
exchange energy-momentum with gravity — locally — yet there is no well 
denned energy-momentum density for gravity itself. This inescapable 
conclusion can be understood as a consequence of the equivalence principle 
(for a discussion see [TB], Section 20.4). 



2. Quasi-local energy- momentum 

The standard approaches aimed at identifying an energy-momentum den- 
sity for gravitating systems always led to various non-covariant, refer- 
ence frame dependent, energy- momentum complexes (such expressions are 
generally referred to as pseudotensors). There are two types of ambigu- 
ity. First, there was no unique expression, but rather many that were 
found by various investigators — including Einstein ES], Papapetrou [33], 
Landau-Lifshitz [IS], Bergmann-Thompson [3J, M0ller [T5], Goldberg [TT] , 
and Weinberg [35] — so which expression should be used? And second — 
in view of the fact all of these expressions are inherently reference frame 
dependent — for a chosen expression which reference frame should be used 
to give the proper physical energy-momentum localization? 

The more modern idea is quasi-local, i.e., energy-momentum should be 
associated not with a local density but rather with a closed 2-surface; for 
a comprehensive review see j2Ej. 

One approach to energy-momentum is via the Hamiltonian (the 
generator of time evolution) . It turns out that this actually includes all the 
classical pseudotensors as special cases, while taming their ambiguities — 
providing clear physical/geometric meaning [il 120). 

3. The covariant Hamiltonian formulation results 

We have developed a covariant Hamiltonian formalism that is applicable 
to a large class of geometric gravity theories [7J El EJ El [201 El EI] ■ For 
such theories the Hamiltonian 3-form H(N) is both a conserved Noether 
current, 

(M{N) oc field eqns ~ 0, (1) 

as well as the generator of the evolution of a spatial region along a space- 
time displacement vector field. It has the general form 

H(N) = N^Hfi + dB{N) , (2) 

where N^H^ — which generates the evolution equations — is itself propor- 
tional to certain field equations (initial value constraints) and thus vanishes 
"on shell" . Consequently the value of the Hamiltonian is determined by 
the total differential (boundary) term: 

E(N,T,):= [ H(N)=<f B(N) . (3) 

Thus, the value is quasi-local. From this boundary term, with suitable 
choices of the vector field on the boundary, one can determine the quasi- 
local energy-momentum and angular momentum/center-of-mass. 



It should be noted that the boundary 2-form B(N) can be modified — by 
hand — in any way without destroying the conservation property. (This is 
a particular case of the usual Nocther conserved current ambiguity.) With 
this freedom one can arrange for almost any conserved quasi-local values. 
Fortunately the Hamiltonian's role in generating evolution equations tames 
that freedom. 

4. Boundary variation principle, reference values 

One must look to the boundary term in the variation of the Hamiltonian 
(see [141 1241 IT5] ) . Requiring it to vanish yields the boundary conditions. 
The Hamiltonian is functionally differ entiable on the phase space of 
fields satisfying these boundary conditions. Modifying the boundary term 
changes the boundary conditions. (The different classical pseudotensors 
are each associated with a specific "superpotential" which can serve as the 
Hamiltonian boundary term, thus they correspond to Hamiltonians with 
different boundary conditions [3].) 

In order to accommodate suitable boundary conditions one must, in 
general, also introduce certain reference values which represent the ground 
state of the field — the "vacuum" (or background field) values. To this end 
for any quantity a we let Acv :— a — a where a is the reference value. 

5. Preferred boundary term for GR 

Some time ago we identified for GR two covariant-symplectic boundary 
terms [7]; one, which was also founcj]] at about the same time by Katz, 
Bicak and Lynden-Bell [HI [12], has become our preferred choiceu 

B(N) = 1 ^(AT a f) Ai N r 1 J + DpN a Ar 1 J), (4) 

This choice corresponds to fixing the orthonormal coframe i? M (equivalcntly 
the metric) on the boundary: 

SH(N)^di N (AT a ASriJ). (5) 

Like other choices, at spatial infinity it gives the ADM [T], MTW [18] , 
Regge-Teitelboim [23], Beig-6 Murchadha g], Szabados energy, 
momentum, angular-momentum, center-of-mass. 

Its special virtues include (i) at null infinity it directly gives the Bondi- 
Trautman energy and the Bondi energy flux [8] , (ii) it is "covariant" , 

1 Via a different route, using a Noether type argument with a global reference. 

2 Here T a p is the connection one-form, n a ^ "' := *(# Q Afl 13 A • • •) and ijv denotes 
the interior product (aka contraction) with the vector field N. 



(iii) it has a positive energy property, (iv) for small spheres it gives a 
positive multiple of the Bel-Robinson tensor, (v) it yields the first law of 
thermodynamics for black holes [3] , (vi) for spherically symmetric solutions 
it has the hoop property [22] . 

6. The reference and the quasi-local quantities 

For all other fields it is appropriate to choose vanishing reference values 
as the reference ground state — the vacuum. But for geometric gravity the 
standard ground state is the non-vanishing Minkowski metric. Thus a 
non-trivial reference is essential. 

Using standard Minkowski coordinates y l , a Killing field of the reference 
has the form N k = Nq + X^iy 1 , where the translation parameters Nq and 
the boost-rotation parameters Ag' = X^ are constants. The 2-surface 
integral of the Hamiltonian boundary term then gives a value of the form 



which yields not only a quasi-local energy-momentum but also a quasi- 
local angular momentum/center-of-mass. The integrals pk(S), Jki{S) in 
the spatial asymptotic limit agree with accepted expressions for these 
quantities ffg] Hfl 13 [27J . 

7. The reference 

For energy-momentum one takes N to be a translation^ Killing field of the 
Minkowski reference. Then the second term in our quasi-local boundary 
expression vanishes. Let us note in passing that holonomically (with 
vanishing reference connection coefficients) the first term in Q reduces 
to Freud's 1939 superpotential [10]. Thus we are in effect here making a 
proposal for good coordinates for the Einstein pseudotensor. 

To construct a reference, choose, in a neighborhood of the desired 
spacelike boundary 2-surface S, four smooth functions y l , i = 0,1,2,3 
with dy° A dy 1 A dy 2 A dy 3 ^ 0; they define a Minkowski reference by 
g = — (dy Q ) 2 + (dy 1 ) 2 + (dy 2 ) 2 + (dy 3 ) 2 . This is equivalent to finding a 
diffcomorphism for a neighborhood of the 2-surface into Minkowski space. 
The reference connection can now be obtained from the pullback of the 
flat Minkowski connection. 

With constant N k our quasi-local expression now takes the form 





B(N) = N' 



k x^ k (T a p-x a ] dy^)M 1 ^ 



(7) 



where dy k = y 



k a dx a has the inverse dx' 



x a k dy k . 



8. Isometric matching of the 2-surface 

The reference metric on the dynamical space has the components 

9iw = Sijy^v 3 v ■ (8) 

Consider the usual embedding restriction: isometric matching of the 2- 
surfacc S. This can be expressed quite simply in terms of quasi-spherical 
foliation adapted coordinates t, r, 9, tp as 

9ab=9ab = 9tjV l AV B = -VaUb + ^HVaVb > ( 9 ) 

where S is given by constant values of t, r, and A, B range over 9, (p. We 
use = to indicate a relation which holds only on the 2-surface S. 

From a classic closed 2-surface into R 3 embedding theorem — as long as 
one restricts S and y°{x fJ ') so that on S 

9ab ■= 9ab + y A y% (10) 

is convex — one has a unique embedding. Wang and Yau have discussed in 
detail this type of embedding of a 2-surface into Minkowski controlled by 
one function in their recent quasi- local work [30| 131] . 

9. Complete 4D isometric matching 

Our "new" proposal is: complete 4D isometric matching on S. (We remark 
that this was already suggested by Szabados back in 200d, and he has 
since extensively explored this idea 26 in unpublished work.) 
Complete 4D isometric matching imposes 10 constraints, 

giiu\s=9nv\s=9ijy t tiV 3 u\s , 

on the 16 y^(to, ?"o, 0, <f)- On the 2-surface S these 16 quantities are 
actually determined by 12 independent embedding functions: y l ,y l t,y l r 
(since from y % on S one can get y l 9,y t v >)- There remain 2 = 12 — 10 
degrees of freedom in choosing the reference. 

One could as an alternative use orthonormal frames. Then the 4D 
isometric matching can be represented by -j? Q =z9 Q . But the reference 
coframe has the form i9 Q = dy a . Thus one should Lorentz transform the 
coframe d a to match dy a on the 2-surface S. This leads to an mtegrability 
condition: the 2-forms dd a should vanish when restricted to the 2-surface: 

^ a |s=0, (11) 



3 At a workshop in Hsinchu, Taiwan. 



This is 4 conditions restricting the 6 parameter local Lorentz gauge 
freedom. Which again reveals that after 4D isometric matching there 
remains 2 = 6 — 4 degrees of freedom in choosing our reference. By the 
way, this orthonormal frame formulation shows that our procedure can 
alternatively be viewed as finding a good frame for the "teleparallel gauge 
current" [5]. 

10. The best matched reference geometry 

There are 12 embedding variables subject to 10 4D isometric matching 
conditions, or equivalently, 6 local Lorentz gauge parameters subject to 
4 frame embedding conditions. To fix the remaining 2, one can regard 
the quasi-local value as a measure of the difference between the dynamical 
and the reference boundary values. So we propose taking the optimal 
"best matched" embedding as the one which gives the extreme value to 
the associated invariant mass m 2 = —piPjg 1 ^ . This is reasonable, as one 
expects the quasi-local energy to be non-negative and to vanish only for 
Minkowski space. 

More precisely, we note two different situations: 
I: Given a 2-surface S find the critical points of m 2 . This should determine 
the reference up to Poincare transformations. 

II: Given a 2-surface S and a vector field N, then one can look to the 
choices of the embedding variables that are a critical point of E(N,S). 
(Afterward one could extremize over the choice of N.) 

Based on some physical and practical computational arguments it seems 
reasonable to expect a unique solution in general. 

For spherically symmetric systems (both static and dynamic), using 
this and some other related strategies we have found reasonable quasi- 
local energy results [33J [TH] |3H |3S] . 
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